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I. INTRODUCTION 



The experimental realization of Bose-Einstein condensates (BECs) in dilute alkali-metal 
atomic vaporsii^ has sparked a large mathematical and physical interest in the study of 
dynamics of condensates, such as solitons^"-, chaos^*^, stability and instabilityii"— , periodic 
and quasi-periodic behaviors^^ii^. 

At ultra-low temperatures, based on mean-field approximation and quasi-one-dimensional 
[quasi-(lD)] regime, the time-dependent condensate wave function ("order parameter") 
iplxyt) is governed by the following cubic nonlinear schrodinger equation (NLS)^^"— 

which is also known as the Gross-Pitaevskii (GP) equation. Here, {ipl"^ is the number density, 
Vo{x) is an external potential, g = [47r/i^a/m][l + C(C^)], and ( = \/[ip^\a\^ is the dilute 
gas parameter. The s-wave scattering length a is determined by the atomic species of the 
condensate. Interactions between atoms are repulsive when a > and attractive when a < 0. 
In coUisionally inhomogeneous BECs, the scattering length is subjected to a spatial periodic 
variation: a{x) = a{x + Lq) for some period Lq leading to nonlinear potentials, which has 
been realized experimentally^"^ and studied theoretically^*^. Also, we can refer to a recent 
comprehensive review^^. 

Spatially periodic potentials V{x), created as optical lattices (OLs), which arise as in- 
terference patterns produced by coherent counterpropagating laser beams illuminating the 
condensate, are of interest in the context of BECs and have been employed in both experi- 
mental and theoretical studies^S"— . 

In this study, we investigate spatially extended solutions of BECs in periodic OLs. We 
apply a coherent structure ansatz to (11. ip . yielding a parametrically forced Duffing equation 
with singularities, which describes the spatial evolution of the field. We employ the averaging 
method to study the periodic orbits (MAWs) including their hyperbolic spatial structures, 
and illustrate their dynamical behaviors with numerical simulation of the GP equation. 

Compared with the previous work— i^"— , the phases of MAWs considered in this paper 
are nontrivial, not constant, which continuously depend on the spatial variation. Nontrivial 
phase solutions have more complicated dynamics^i^ and imply nonzero current of the matter 
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- it is proportional to E?{x)9{x) = c, for amplitude R{x) of MAW and nonzero constant 
c -along X-axis, and hence seem to have no direct relation to present experimental setting 
for BEC9^'^(remember that the condensate in this paper is confined to be a parabolic 
trap). Of course, dealing with the MAWs with nonzero phases will take more difficulties 
than before. One reason is that the forced Duffing oscillator derived from GP equation has 
singularities at the origin and the directly application of the usual perturbation theory is 
unavailable because of the loss of smoothness and the existence of strong nonlinear term 
(singular term). 

The averaging method^""— at its heart is a transformation procedure leading to a sys- 
tematic perturbation expansion and completes with error bounds on the difference between 
exact and approximation solutions. Also, it is a important tool for proving properties of 
the exact problem based on properties of the approximation problent^. For example, the 
existence of periodic orbits can be proved using averaging together with the implicit function 
theorem, and the existence of invariant tori can be proved using averaging together with the 
Moser twist theorem. However, applying the averaging method for singular system, the first 
problem we must deal with is how to put the system into the standard form of averaging. 

The paper is organized as follows. In Section [Til we introduce modulated amplitude 
waves involving periodic and quasi-periodic, and in Section IIIH we apply a transformation 
to transform the GP equation to a standard form of averaging. An abstract result upon the 
dynamics to BECs is obtained in Section IIV[ and in Section |V] we analyze and demonstrate 
some of the spatial dynamical features of BECs with a positive chemical potential. Finally, 
we summarize our results in Section I VII 

II. COHERENT STRUCTURE AND MODULATED AMPLITUDE WAVE 

We consider uniformly propagating coherent structures with the ansatz 

^{t,x) = R{x) exp{i[e{x) - nt]), (2.1) 

where R{x) G M gives the amplitude dynamics of the condensate wave function, 6{x) de- 
termines the phase dynamics, and the "chemical potential" /i, defined as the energy which 
takes to add one more particle to the system, is proportional to the number of atoms 

3 



trapped in the condensate. When the (temporally periodic) coherent structure (12. ip is also 
spatially periodic, it is called a modulated amplitude wave (MAW)^'^. Similarly, a solu- 
tion of the equation (II. ip with the (temporally periodic) coherent structure (12. ip is called a 
quasi-periodic modulated amplitude wave (QMAW) if it is also spatially quasi-periodic. 

Inserting (12. ip into (II. ip . we obtain the following two couple nonlinear ordinary differen- 
tial equations 

R" + 5R - + eaR^ + eV{x)R = 0, (2.2) 

e" + 2Q'R'/R = ^ Q'{x) = (2.3) 



where 



5 := ea := — eV{x) := -^Vo{x) 



and the integration constant c, determined by the velocity and number density, plays the 
role of "angular momentum"—. 

Inspecting (12. 2p we know that in case of c = 0, i.e., the phase of the condensate wave 
function is trivial, it is the parametrically driven Duffing equation with the time variable 
replaced by the spatial coordinate, and MAWs (standing waves) in this system with V{x) = 
Vo cos Kx or V{x) = V\ cos k^x + Vi cos k^x have been widely studied^^'^'^. 

In general, c 7^ 0, the system (12. 2p becomes more complicated and the phase is no longer 
constant^. Even the amplitude R{x), a solution of (12. 2p . is L-periodic, the corresponding 
condensate wave function ■iIj{x, t) may be not periodic with respect to the spatial variable x. 
In fact, 

ip(t, x) = R{x)expi[Q{x) — fit] 

= R{x) ( cos[6(x) + ux — fit] + i sin[6(x) + ux — /it]) , 



1 r°^^ c 



where 

and 0(x) = G(x) — i/ is a L-periodic function with zero mean value. If 2'it/u and L are 
rationally related, then ip{x, t) is a MAW; if 27r/i^ and L are rationally irrelevant, then ip{x, t) 
is not periodic but quasi-periodic, which is corresponding to a QMAW with the frequency 
u = {2-71 /u, L). 



There also exists an interesting and surprising result. Note that i/ = when c = and 
u ^ when c 7^ 0. If we vary c on the interval (—00, +00), by continuous dependence of 
solutions with respect to the parameters, u can continuously take the value on some interval, 
which implies that (11. ip has infinitely many (positive measure set) MAWs and QMAWs by 
adjusting the integration constant c. Thus, all these MAWs and QMAWs form a foliation. 

In this paper, we consider the case 6 > corresponding to a positive chemical potential. 
Also, in order that the mathematical results obtained in this paper do apply to more general 
periodic functions, we assume that the external potential V{x) is an analytic and L-periodic 
function (OLs). Note that (12. 2p defines on two half-planes, and we only consider the case of 
the right half-plane since there are no distinct technicalities. 



III. TRANSFORMATION TO STANDARD FORM OF AVERAGING 



Rewrite equation (12.20 in the planar equivalent form 



R' = S 

S' = -6R + --^ - eaR^ - eV{x)R. 
R^ 



(3.1) 



Generally, averaging method involves two steps: transforming to standard form; solving the 
averaging equation. In order to proceed we need to transform (13. ip to a standard form for 
the method of averaging. So we have the following result. 

Lemma 3.1. Under the transformation \1/ : T x (y^, +00) (0, +00) x M defined by 



R = p\ cos'^{\/6x 



+ ^sin^(v^x + ^) 
op* 

cos(-\Ax + 6) sm{\/5x + 



' cos^ 



sin^^y/Sx + 
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system (13. ip changes into a new system 



a 



9' = s 



1 1 g2 

-(1 + j^) sm2{VSx + 9) + -(1 - — ) sin4(v^x + 9) 
+ -^V{x)sm2{VSx + 9)\ 

(3 + cos i{VSx + #)) + 3°^ty/_'ll) "^"^ 2(^^ + ") 



+ ^,V(.) 1 1 + 



wt/i the new coordinates {9, p) in the half-plane T x (-^ y, +oo). 

Proof. First, it is easy to verify that, for each p G (y +oo) and ^ G T, 



(3.2) 



(/9(x; 6*, p) = p\ cos'^{V6x + 9) + -^—r sm'^{V6x + 9) 
V '^P 

is a periodic solution with the same period r = 2n / \^ of the unperturbed equation 

R" + 6R-4^ = 
or the equivalent Hamiltonian system 



with the initial value 



R' = S 



<^(--^;^,p)=p, <^'(--^;e,p) = 0. 



(3.3) 



That is to say, (13. 3p is an isochronous system. The total energy of system (13.30 is given by 

\s' + V{R) = V{pl 



where 



V{R) 



5R^ c2 



2 2i?2' 

and as we know, all periodic solutions lie on curves of constant energy. We will use these 
facts below. 
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Now using the variation of constants, the functions p(-), 6{-) can be defined such that 

dip 

Rix) = (p(x + 9(x), p(x)), S(x) = —(x + 6(x),p(x)). 

ox 

From the conservation of the Hamiltonian, it follows that 

^(^)^ + ^<^) = 

so the differentiating with respect to x along with system 03.11] yields 
,,,^^dp d /I 2, 



e [aR^{x, e, p) + V{x)R{x, 9, p)] S{x, 6, p). 



Then, we have 



So, using the definition of \1/ together with equation (13. 5p . we have the first desired expression 
of (O- 

Using the formula for B! given in system (13. ip and the definition of p(-) and 6'(-), we 
obtain that 

dip f ^ , '^^ ^'■P 



dx \ dx J dp dx dx 
Together with equation (13. 5p . after some simple algebraic manipulations, it follows that 

£ = [^R'i^^ p) + vi^M^^ p)] 

Finally, using the definition of ^ together with equation (13. 5p . we have the second desired 
expression of (13. 2p . □ 
In general, finding a explicit expression of the transformation is not easy. The transfor- 
mation \1/ given in Lemma 13.11 has a more delicate information, such as protecting two- form 

di? A d5 = y/5{p - -^)dp A d^, 
bp'^' 

which is not used in this paper. However, we believe that it will be helpful for further study. 



IV. AN ABSTRACT RESULT OF AVERAGING TO BEGS 



In this section, we will give an abstract result to BECs by the method of averaging. We 
also assume that the period L of the external potential V{x) satisfies that 

- = -eQ, (g,p) = l, p,geZ+, 
r p 

where r is the least period of the unperturbed system. The basic idea that leads to the 
application of the method of averaging arises from an inspection of system f l3.2p . The 
derivatives with respect to the spatial variable x are all proportional to e. Hence, if e is 
small, the variables would be expected to remain near their constant unperturbed values 
over a long spatial scale. 

A good approximation of f l3.2p up to the spatial domains of order 1/e is given by the 
averaged system 



where 



p' = ep^{d) 



^{9) = [ V{s) sin 2(755 + 9)ds, 

2v 6L Jo 

3a 

ao 



86^ L 

ai = i / V{s)ds, L = qpmm{L,T}. 
2dL Jo 

We have the following theorem. 

Theorem 4.1. There exists a c'^,r > 2, change of variables 

p = p + ewi{6,p,x,e), 6 = 6 + ew2{0,p,x,e) 

with wi,W2 L -periodic functions ofx, transforming ( 13. 2p into 
p' = ep<l>{9)+e^gi{9, p,x,e) 
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(4.1) 



'withgi,g2 L-periodic functions ofx. Moreover, 

(i) If {9e{x), pe{x)) and {9q{x) , po{x)) are solutions of the original system (13.21) and the 
averaged system (14.11) respectively, with the initial values such that 

\ pM - Pom + \ 0e{0) - 9om = 0{e), 

then 

I p,(x) - poix)\ + I 9,ix) - 9oix)\ = Oie), 

for the spatial domains x of order 1/e. 

(ii) // there exist two constants po ^ i\f^y +00), 6'o G M such that 

$(^0) = 0, (4.2) 

^{6-pt + c') + a,+ nOo) = 0, (4.3) 

<t>\9o) ^ 0, (4.4) 

2ao{6po - ^) - ^^0§^,nOo) + 0, (4.5) 

i.e., {9o,Po) is an equilibrium point of (14.11) such that the corresponding Jacobian matrix has 
no eigenvalue equaling to zero, then i\3.2h admits a L-periodic solution {9e{x), ps{x)) such 
that 

\pe{x)- po\ + \9,{x)-9o\ = 0{e), 

for sufficiently small e; if, in addition, 

.2 4v^c2'^3 



2aoi5po - ^Wi9o) - ,,_4 ° > 0, (4.6) 

Po i^Po - c^r 



then the L-periodic solution {ps{x), 9^{x)) is hyperbolic and instable with respect to the spatial 
variable x. 

Proof. Conditions (l4.2p -f Hy6|l imply that {9q, po) is a instable and hyperbolic fixed point 
of system (14. ip . The proof of this theorem follows directly front^ (Theorem 3.2.3) or— 
(Theorem 4.1.1). □ 

Remark 4.1. The instability in Theorem 4-i is only relevant to the amplitude equation 
(2.2), which is some artificial "instability" in terms of the evolution in x. There is a set of 
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methods for the study of modulational instability in time t, e.g., se^~—. Recently, based on 
spectrum theory and Hamiltonian floquet theory, the method of studying modulational ( tem- 
poral) instability of standing waves (with trivial phases) has been developed for NLS equation 
with constant nonlinearity coefficients and periodic potentials by Bronski and RaptM-, later 
applied by Porter et alM. However, here we can not provide any information upon it with 
our methods. 



V. EQUILIBRIUMS AND THE AVERAGED EQUATION 

In this section, we will analyze and demonstrate some of the spatial dynamical features 
of BECs with a positive chemical potential. We also assume that V{x) is an analytic and 
L-periodic function with the least positive period L = r = ir/y/S, i.e., p = q = 1. 

First, expanding V in a Fourier series, we have 

oo 

V{x) =bo + ^[ofc sm{2V6kx) + bk cos{2^/Skx)], (5.1) 

k=l 

where all coefficients are real. Let us substitute the expansion for V{x) given by (15. ip into 
(14. ip . and by an easy (perhaps lengthy) computation, we obtain the averaging system 




2 + bj sm{2e 



4v^ , , 

(5 2) 

, _ 3ea{6p^ + c^) ebo e 6p^ + ' ^ ' ^ 



where 



(j) = arctan ^ {(f) = Sign(ai) ■ 7r/2, if bi =0). 
bi 

If a < 0, corresponding to the repulsive nonlinearity, we take 6q = kn — 0/2. Let us 
define the function / : ( +oo) — )• M by 

3a{ 6p^ + c^) ebo e_8fW^ 
^ 25 ^ 45 5p4_c; 



It is easy to verify that 

/(p) — )■ — oo, as p — 7- +oo; /(p) — +oo, as p — )■ \J (? jb 



■+ 
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Using the mean- value theorem, there is at least one root po of / on the interval ( y +00). 
If a > 0, corresponding to the attractive case, we take 6q = kir + {n — 0)/2. Similarly, we 
can induce that the function / also has at least one root po on the interval +00). 
Thus, in any case, the averaged system (15.21) has at least one equilibrium (^0;Po)- 

After a simple computation, the Jacobian matrix of the averaged system fl5.2p at the 
equilibrium (6'o, po) is given by 



M 



( 



\ 



ep^)\/a\ + h\ cos(2^o + 0) \ 







(5.4) 



Spl{6pl + c^){6pl-c') 

Notice that if M has no zero eigenvalue, then the equilibrium {6q, po) can be continuable. If 
M has no imaginary eigenvalue, the equilibrium {6q, po) is hyperbolic since the two eigenval- 
ues of M have opposite signs. So it follows that the original system (13. 2p has a hyperbolic 
L-periodic solution (^^(x), Pe(x)) with respect to the spatial variable x near (6'o,Po) such 
that {9^{x),ps{x)) (6*0, Po), as e 0. 

To demonstrate the process of averaging to BECs, a specific example of numerical compu- 
tation is given in the following. We take the integration constant c = 1 and the parameters 
6 = l,a = 0.3,60 = — 12.4,61 = — 20.19, ai = 0. By a numerical computation via MATH- 
METICA, we obtain equilibriums 

(fc7r + |,3.00), (A;7r + |,2.00), (A;7r, 10.00), A; = 0, ±1, ±2, • • • 



for the averaged system (15. 2p . 

As we know, {Oq,pq) = (A;7r -|- 7r/2, 3.00) and (6*1, pi) = {kn, 10.00) are hyperbolic with the 
eigenvalue of linearization 

Ao = ±0.95e, Ai = ±15.33e. 

By the averaging theorem, the equilibriums can persist as the periodic orbits for the original 
system (13.11) : in addition, these periodic orbits are also hyperbolic with respect to spatial 
variable x. 

Returning to (13. ip . consider its unperturbed system, and using Lemma [3TT| the periodic 
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FIG. 1. A plot of a MAW ipQ{x,t) corresponding to the equilibrium 3) for unperturbed system 
of (jl.ip on the right-plane. The original system (jl.ip has a MAW or QMAW ip£{x,t) such that 
I |'(/'o(2^) i)| ~ |V'e(^5^)l| = C'(e) for all time t and the spatial domains of order 1/e. The parameters 
we take as follows: c = 1, 6 = 1, a = 0.30, 6o = —12.40, bi = —20.20. (a) Spatial amplitude R{x) 
plot; (b) Nontrivial spatial phase Q{x) plot; (c) A plot of space-time Re[ijj{t,x)] ; (d) The density 
plot of Re[ip{t, x)]. 

orbit corresponding to the equilibrium (6'o,Po) = {kn + 7r/2,3.00) is given by 



(po{x) = SWsin X + — cos^ x. 



Following from (12. 3p . we have the angle function 



eo(a;) 



ds 
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FIG. 2. A plot of a MAW ipo{x, t) corresponding to the equilibrium (vr, 10) for unperturbed system 
of (jl.ip on the right-plane. The original system (jl.ip has a MAW or QMAW ipe{x,t) such that 
I |'0o(3;ii)| — |V'£(^i^)l| = ^(^) foi' time t and the spatial domains of order 1/e. The parameters 
we take as follows: c = 1, 5 = 1, a = 0.30, 6o = —12.40, &i = —20.20. (a) Spatial amplitude R{x) 
plot; (b) Nontrivial spatial phase Q{x) plot; (c) A plot of space-time Re[^/'(t,x)] ; (d) The density 
plot of Re[i/'(t, x)]. 



with its mean value 



eo 



9(sin^ X + ^ cos^ x) 



ds = 1. 



We conclude that ipoit, x) = Ro{x) exp{i[Qo{x) — fit]) is a MAW of unperturbed system (11. II) . 
see FlCm The averaging theorem implies that the MAW ilJo{t,x) can be continuable, i.e., 
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there exists a MAW or QMAW tjj£{t,x) for system (11. ip such that 



\\ije{t,x)\-\Mi,x)\\=0{e), 

for sufficiently small e. We can have a similar analysis for the equilibrium (6i, pi) = 
(A;7r, 10.00), see FIG. El 

The equilibrium {62, P2) = (A;7r + 7r/2, 2.00) with the linearized eigenvalue A2 = =F6.44ei 
also persist as periodic orbits for system (13. ip . Since the equilibrium is not hyperbolic, one 
cannot depict the spatial dynamics of the corresponding continuable periodic orbit. This 
question is left open for further study. 

VI. DISCUSSION AND CONCLUSION 

We have given an abstract result to BECs, see Theorem 14.11 Using averaging method, we 
determine the location where the MAWs or QMAWs exist. Comparing the previous work, 
we do not restrict our discussion near the origin since the equation we have dealt with has 
some singularities. On the other hand, the MAWs or QMAWs studied in this paper, which 
have nontrivial phases, form a foliation. This is a new result in this context. 

However, we can not determine the spatial dynamics of some periodic orbits corresponding 
to MAWs or QMAWs since the equilibriums are not hyperbolic. Maybe second or even 
higher-order averaging is required, and this question is left for our further study. In addition, 
we can not provide any information upon modulational instability in time t of MAWs or 
QMAWs. One reason is that MAWs or QMAWs obtained in this paper are not standing 
waves, but the waves with nontrivial phases, which does not allow us to apply the method 
developed by Bronski and Rapti^ directly. This problem may be not easy but a good topic 
for further study. 

From the view point of a physical application, it might be reasonable to use the averaging 
principle to replace a mathematical model by the corresponding averaged system, to use the 
averaged system to make a prediction, and then to test the prediction against the results of 
a physical experiment. The study of this paper exactly gives a frame for application of the 
averaging method to BECs. 
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